Emergence of Complex Spatio- Temporal Behavior in Nonlinear Field Theories 
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We investigate the emergence of time-dependent nonperturbative configurations during the evolu- 
tion of nonlinear scalar field models with symmetric and asymmetric double-well potentials. Com- 
plex spatio-temporal behavior emerges as the system seeks to establish equipartition after a fast 
quench. We show that fast quenches may dramatically modify the decay rate of metastable states 
in first order phase transitions. We briefiy suggest possible applications in condensed matter systems 
and early universe cosmology. 
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INTRODUCTION 

The emergence of complex patterns is one of the most 
distinct signatures of nonlinear interactions in natural 
systems. Since Einstein's pioneering work on Brownian 
motion [l|, it became clear that much can be accom- 
phshed by modeling the interactions of a system with its 
environment through the action of random and viscous 
forces. During most of the twentieth century, studies 
were mainly restricted to investigating the motion of a 
point particle in nonlinear potentials ^. With the ad- 
vent of fast computers, modeling of stochastic evolution 
added spatial dimensions, allowing for the quantitative 
study of spatio-temporal complex behavior. Up to about 
ten years ago, most of the work concentrated in hydro- 
dynamical and soft condensed-matter systems [2|. Re- 
cently, developments in high energy physics and cosmol- 
ogy have opened the interesting possibility that complex 
spatio-temporal behavior may also play a role in relativis- 
tic field theories, in particular during the early stages of 
cosmological evolution |j| and may even be observed in 
high-energy colliders |3|- 

Here we will briefly review some of the work done dur- 
ing the past few years which focused on understanding 
the effects of fast quenches on nonlinear scalar field the- 
ories. The quenches model both temperature quenches 
in the context of fast cosmological expansion (in partic- 
ular at scales close to the GUT scale ~ 10^^ GeV) or the 
cooling of fireballs during high energy collisions such as 
those currently being investigated at RHIC and soon at 
LHC. The quench may also represent a pressure quench, 
common in condensed matter physics or, more generally, 
the appearance of a low-energy effective interaction that 



modifies the effective potential of the long-wavelength 
modes of the field or order parameter describing the sys- 
tem's evolution. 

THE MODEL 

Consider a (2+l)-dimensional real scalar field (or 
scalar order parameter) (/'(x, t) evolving under the influ- 
ence of a potential V{(j)). The continuum Hamiltonian is 
conserved and the total energy of a given field configura- 
tion (?!)(x, t) is, 
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where V(4i) — ^cp^ — ^cj)^ + ^<j)* is the potential energy 
density. The parameters m, a, and A are positive definite 
and temperature independent. It is helpful to introduce 
the dimensionlcss variables cj)' = (pyX/m, x' — xm, t' — 
tm, and a' = a/{m^/X) (Wc will henceforth drop the 
primes). Prior to the quench, a = and the potential is 
an anharmonic single well symmetric about (f) — 0. The 
fleld is in thermal equilibrium with a temperature T. At 
the temperatures considered, the fluctuations of the fleld 
are well approximated by a Gaussian distribution, with 
((/)^) — aT (a — 0.51 and can be computed numerically). 
As such, within the context of the Hartree approximation 
|6j, the momentum and field modes in fc-space can be 
obtained from a harmonic effective potential, and satisfy 
(|7f(fc)|2) = T and (|0(k)|2) = ^^, respectively. The 

Hartree mass rnjj = 1 + ^{4>^) depends on the magnitude 
of the fluctuations (and thus T). Within the Hartree 
approximation we can write the effective potential as 
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FIG. 1: Schematic picture showing change in potential V{<l)) 
from single-well to symmetric double-weU after quench. 



Hereafter we will refer to a particular system by its initial 
temperature. All results are ensemble averages over 100 
simulations. 

If a ^ 0, the Z2 symmetry is explicitly broken. When 
a = 1.5 = ttc, the potential is a symmetric double-well 
(SDW), with two degenerate minima. This is the first 
case we consider. 



QUENCHING INTO SYMMETRIC DOUBLE 

WELLS: EMERGENCE OF SPATIO-TEMPORAL 

ORDER 

At a — Uc — 1.5, the quench amounts to switching 
from a single to a double well with the field initially lo- 
calized at = 0. In Fig. [Jwe indicate this schematically. 

As shown in Ref. 10], the quench sets oscillations 
in the field's zero mode, (/)avc = y J'f'dV, where V is 
the volume (or area in 2d). The amplitude of these os- 
cillations is controlled by the temperature of the initial 
Gaussian distribution, as explained above. In fact, tem- 
perature here is simply a convenient way to set an initial 
Gaussian distribution in momentum space. We did this 
using a Langevin equation with white noise. One could 
state that the field is at T = but initially set with a 
Gaussian distribution in momentum space with a certain 
width. This width is a measure of the initial "tempera- 
ture" of the system. 

At early times small fluctuations satisfy a Mathieu 
equation in fc-space 

Scf>=-[k^ + V:if[^,^,{t)]]Scl,, (3) 

and, depending on the wave number and parametric os- 
cillations of (f>s^vc{t), can undergo exponential amplifica- 
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FIG. 2: Number of oscillons nucleated between t and t + St a,t 
T = 0.22 and St = 1. The global emergence is evident early 
in the simulations. Inset: probability distribution of radii and 
periods of oscillations of individual oscillons. 
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For T < 0.13, no modes are ever amplified. 



As the temperature is increased, so is the amplitude and 
period of oscillation in i^ave, gradually causing the band 
< fc < 0.48 to resonate and grow. Furthermore, for 
large enough temperatures (T > 0.13) large-amplitude 
fluctuations about the zero mode probe into unstable re- 
gions where V^g < 0, which also promote their growth. 
Note that this is very distinct from spinodal decomposi- 
tion, where competing domains of the two phases coarsen 
[3. Instead, for the values of T and a considered, 0avc 
continues to oscillate about the (^ = minimum. 

As a result of the energy transfer modeled by para- 
metric amplification, oscillons are nucleated initially in 
phase. But what are oscillons? They are the higher di- 
mensional equivalent of kink-antikink breathers, familiar 
of Id nonlinear dynamics [Sj. Extensive work has been 
done on oscillons and their properties and the reader can 
consult the relevant literature listed in Ref. [9j. Here, it 
is enough to mention that oscillons are long-lived, time- 
dependent, localized field configurations which express 
local ordering of momentum modes. What was also ob- 
served in Ref. [li| is that after the quench oscillons 
emerge in synchrony, exhibiting both spatial and time 
ordering. In Figure Owe illustrate this phenomenon. 

Finally, we introduce a measure of the partitioning of 
the kinetic energy H(i), which we use to describe the 
nonequilibrium evolution of the system: 
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where p(k, t) = K(k, t)/ J (fkK(k, t), and K{k, t) is the 
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FIG. 3: The change of n(i) from the initial state for closed 
systems at T = 0.22. The exponential approach to equilib- 
rium is clear at late times. The inset illustrates the role of 
oscillons as a bottleneck to equipartition. 



kinetic energy of the fc-th mode. Hit) attains its maxi- 
mum (n,„ax = In(A^) on a lattice with N degrees of free- 
dom) when equipartition is satisfied. This occurs both 
at the initial thermalization (t — 0) and final equilib- 
rium states, since in this case all modes carry the same 
fractional kinetic energy. In Fig. |3|we show the change 
of n(t) from the initial state, n(i = 0) — n(t), for the 
closed system at T = 0.22. At late times {t > 150), we 
have found that the system equilibrates exponentially in 
a time-scale r ~ 500. At early times, the localization 
of energy at lower k-modes, corresponding to the global 
emergence of oscillons, prolongs this approach to equipar- 
tition. The inset of Fig. |2| shows the large variations in 
n(t) (dotted line) that arise due to the synchronous oscil- 
lations in the kinetic energy of these configurations. Also 
shown (solid line) is the average between successive peaks 
of n(i), with a plateau at approximately 20 < t < 70 
that coincides with the maximum oscillon presence in 
the system. Thus, oscillon configurations serve as early 
bottlenecks to equipartition, temporarily suppressing the 
diffusion of energy from low (0 < |k| < 0.8) to higher 
modes. 



QUENCHING INTO ASYMMETRIC DOUBLE 
WELLS: RESONANT NUCLEATION 

For a > Uc — 1-5 the potential is asymmetric with the 
minimum at = becoming metastable. We proceed 
as before by quenching the system from a single well, as 
illustrated in Fig. ^ 

As before, oscillons will once again be nucleated. How- 
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FIG. 4: Schematic picture showing change in potential V^(0) 
from single-well to asymmetric double-well after quench. 



ever, the situation now changes dramatically. Due to the 
asymmetry, the system will decay into the global mini- 
mum at 0+ . We have observed that this decay may occur 
in three possible ways depending on the initial temper- 
ature T and the value of a [llj : i) the transition to the 
global minimum happens very fast in what can be called 
a "cross-over" transition; ii) the transition occurs as a 
single oscillon becomes unstable and grows into a critical 
bubble. As is well know from the theory of first order 
phase transitions |3| , once a critical nucleus forms it will 
grow to complete the transition; iii) two or more oscillons 
coalesce to become a critical bubble that then grows to 
complete the transition. 

In order to simplify the analysis, we fixed the temper- 
ature to be r < 0.22. From the Hartree potential of 
Eq. 121 one can see that for large temperatures the po- 
tential becomes a single well again. For T < 0.13 no 
oscillons are nucleated after the quench. In this case, 
we expect that the usual metastable decay rate based 
on the theory of homogeneous nucleation (HN)will ap- 
ply, becoming more accurate for smaller T [3, Il2| | . The 
decay rate per unit volume obtained from HN theory 
is controlled by the Arrhenius exponential suppression, 
r(r,a) ~ T('^+i)exp[-^6(T,a)/T], where Et is the en- 
ergy of the critical bubble or nucleus and d is the number 
of spatial dimensions. [We use units where c = ks = h = 
1.] The typical time-scale for the decay in a volume V is 
then, THN ^ (VT)-^ ^ T-^exp[Eb{T,a)/T]. 

In Fig. [Slwe show the evolution of the order parameter 
4'aNc{i) as a function of time for several values of asymme- 
try, 1.518 <a< 1.746, for T = 0.22. Not surprisingly, as 
a — > ttc = 1.5, the field remains longer in the metastable 
state, since the nucleation energy barrier E^ ^ oo at 
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FIG. 5: The evolution of the order parameter <jf>avc(i) at 
T = 0.22 for several values of the asymmetry. From left to 
right, a = 1.746, 1.56, 1.542, 1.53, 1.524, 1.521, 1.518. The in- 
set shows Vcff for the same values. 
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FIG. 6: Decay time-scale trn as a function of critical nucle- 
ation effective free-energy barrier Ef,/T at T = 0.18, 0.20, and 
T = 0.22. The best fits (dashed lines) are power-laws with 
exponents B ~ 3.762, 3.074, and 2.637, respectively. 



ttc- However, a quick glance at the time axis shows the 
fast decay time-scale, of order 10^~^. For comparison, 
for 1.518 < a < 1.56, HN would predict nucleation time- 
scales of order - lO^^ > thn ~ exp[Eb/T] > lO^^ (i^ ^[_ 
mensionless units). [The related nucleation barriers with 
the effective potential are £";,(« = 1.518) = 14.10 and 
Ei,{a = 1.56) = 5.74.] For small asymmetries 0avc(i) dis- 
plays similar oscillatory behavior to the SDW case before 
transitioning to the global minimum. As a is increased 
the number of oscillations decreases. For large asym- 
metries, a > 1.746, the entire field crosses over to the 
global minimum without any nucleation event, resulting 
in oscillations about the global minimum. This is the 
situation described in case i) above. 

In Fig. we show the ensemble-averaged nucleation 
time-scales for resonant nucleation, trn, as a function of 
the nucleation barrier (computed with Eq. [5J), Ei,/T, for 
the temperatures T = 0.18, 0.20, and 0.22. [For tem- 
peratures above T = 0.26 we are in the vicinity of the 
critical point in which no barrier exists.] The nucleation 
time was measured when 0ave crosses the maximum of 
Vcs. The best fit is a power law: 
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with B ^ 3.762 ± 0.016 for T = 0.18, B = 3.074 ± 0.015 
for T = 0.20, and B = 2.637 ± 0.018 for T = 0.22. This 
simple power law holds for the same range of tempera- 
tures where we have observed the synchronous emergence 
of oscillons. It is not surprising that the exponent B in- 
creases with decreasing T, since the synchronous emer- 
gence of oscillons becomes less pronounced and eventu- 
ally vanishes. In these cases we should expect a smooth 
transition into the exponential time-scales of HN. 



We conclude that fast quenching can dramatically af- 
fect the nucleation time-scale of first order phase transi- 
tions. In other words, HN fails for fast quenches. 

Here we propose the mechanism by which this fast de- 
cay occurs: for nearly degenerate potentials, ac < a < 
ai, the critical nucleus has a much larger radius than a 
typical oscillon; it will appear as two or more oscillons 
coalesce. We call this Region I, defined for Rj, > 2i?osc, 
where RoscJ^s the minimum oscillon radius computed 
from Ref. [a]- Figure illustrates this mechanism. Two 
oscillons, labeled A and B, join to become a critical nu- 
cle us. [The interested reader can see simulation movies 
at jhttp://www. dartmouth.edu/~cosmos/oscillons] 

As a is increased further, the radius of the critical nu- 
cleus decreases, approaching that of an oscillon. In this 
case, a single oscillon grows unstable to become the crit- 
ical nucleus promoting the fast decay of the metastable 
state: there is no coalescence. We call this Region II, 
ai < a < au, Rb < 2i?osc- This explains the small num- 
ber of oscillations on (j>ave{t) as a is increased [cf. Fig. 
|3]. To corroborate our argument, in Fig. |Slwe contrast 
the critical nucleation radius with that of oscillons as ob- 
tained in Ref. [5j , for different values of effective energy 
barrier and related values of a at T — 0.22. The crit- 
ical nucleus radius Rb is equal to 2i?osc for a — 1.547. 
This defines the boundary between Regions I and II: for 
a > ai a single oscillon may grow into a critical bubble. 
Finally, for a > an = 1.746 the field crosses over to the 
global minimum without any nucleation event. 

An obvious extension of the present work is the in- 
vestigation of "resonant nucleation" in 3d. Preliminary 
results indicate that the power law behavior persists with 
_B ~ 1.5 for the relevant range of temperatures for oscil- 




FIG. 7: Two oscillons coalesce to form a critical bubble. First 
two frames from top show oscillons A and B. Third and fourth 
frames shows A and B coalescing into a critical bubble. Fi- 
nal frame shows growth of bubble expanding into metastable 
state. 



Ion coalescence. An application in cosmology has recently 
been sug gested in the context of a two-field model of infla- 
tion [13. In general, RN will occur whenever the effective 
potential changes faster than the typical relaxation rate 
of the longest wavelength of the order parameter. These 
results could be extended to systems in the Ising uni- 
versality class, in particular to ferromagnetic transitions, 
where the quench could be induced by a homogeneous 
external magnetic field. These topics are presently under 
investigation. 

One of us (MG) would like to thank the organizers, 
in particular Jean-Michel Alimi for being a wonderfully 
gracious host. 




FIG. 8: Radius of critical bubble (Rt) and twice the minimum 
oscillon radius (2_Rosc) as a function of its energy barrier and 
related values of a at T = 0.22. For a > 1.547 one cannot 
easily distinguish between an oscillon and a critical bubble. 
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